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Abstract We review the mechanisms via which an external magnetic field can affect 
the ground state of cold and dense quark matter. In the absence of a magnetic field, 
at asymptotically high densities, cold quark matter is in the Color-Flavor-Locked 
(CFL) phase of color superconductivity characterized by three scales: the supercon- 
ducting gap, the gluon Meissner mass, and the baryonic chemical potential. When an 
applied magnetic field becomes comparable with each of these scales, new phases 
and/or condensates may emerge. They include the magnetic CFL (MCFL) phase 
that becomes relevant for fields of the order of the gap scale; the paramagnetic CFL, 
important when the field is of the order of the Meissner mass, and a spin-one con- 
densate associated to the magnetic moment of the Cooper pairs, significant at fields 
of the order of the chemical potential. We discuss the equation of state (EoS) of 
MCFL matter for a large range of field values and consider possible applications of 
the magnetic effects on dense quark matter to the astrophysics of compact stars. 



1 Introduction 

Contrary to what our naive intuition might indicate, a magnetic field does not need 
to be of the order of the baryon chemical potential to produce a noticeable effect 
in a color superconductor. As discuss in ll46l [62], a color superconductor can be 
characterized by various scales and different physics can occur at field strengths 
comparable to each of them. Specifically, for the so-called Color-Flavor-Locked 
(CFL) phase, the superconducting gap, the Meissner mass of the charged gluons, 
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and the baryon chemical potential define three scales that determine the values of the 
magnetic field needed to produce different effects. Thus, the presence of sufficiently 
strong fields can modify the properties of the dense-matter phase which in turn 
might lead to observable signatures. 

In this paper we review the status of our current knowledge about the magnetic 
field effects on color superconductivity (CS) at asymptotically high densities and 
discuss possible consequences of these effects for the physics of compact stars. 



2 Magnetic Fields in Compact Stars 

The density of matter in the core of compact stars is expected to exceed that of nu- 
clear matter p nuc = 2.8 x 10 14 g/crrr 1 1 128 1. At such densities, individual nucleons 
overlap substantially. Under such conditions matter might consist of weakly inter- 
acting quarks rather than of hadrons 11371 . Due to the asymptotic freedom mecha- 
nism l36l l76l II 191 one might think that at high baryon density QCD is amenable 
to perturbative techniques [19] [22] [70] ED [88] [S9] [TD5|. However, the ground 
state of the superdense quark system, a Fermi liquid of weakly interacting quarks, 
is unstable with respect to the formation of diquark condensates lfT7ll2Tll69l , a non- 
perturbative phenomenon essentially equivalent to the Cooper instability of BCS 
superconductivity. Given that in QCD one gluon exchange between two quarks is 
attractive in the color-antitriplet channel, at sufficiently high density and sufficiently 
small temperature quarks should condense into Cooper pairs, which are color an- 
titriplets. These color condensates break the SU(3) color gauge symmetry of the 
ground state producing a color superconductor. 

In the late 90's the interest in CS was regained after the finding, based on dif- 
ferent effective theories for low energy QCD [5 122J, that a color-breaking diquark 
condensate of much larger magnitude than originally thought may exist already at 
relatively moderate densities (of the order of a few times the nuclear matter density). 
At densities much higher than the masses of the u, d, and s quarks, one can assume 
the three quarks as massless. In this asymptotic region the most favored state is the 
CFL phase (5|, characterized by a spin-zero diquark condensate antisymmetric in 
both color and flavor. Since the combination of high densities and relatively low 
temperatures can naturally exist in the dense cores of compact stars, it is expected 
that CS could be realized in that astrophysical environment. 

Compact stars, on the other hand, are very magnetized objects. From the mea- 
sured periods and spin down of soft-gamma repeaters (SGR) and anomalous X- 
ray pulsars (AXP), as well as the observed X-ray luminosities of AXR it has been 
found that a certain class of neutron stars named magnetars can have surface mag- 
netic fields as large as 10 14 - 10 16 G (83] [92] [109] [137] [[38) . In addition, since the 
stellar medium has a very high electric conductivity, the magnetic flux should be 
conserved. Hence, it is natural to expect an increase of the magnetic field strength 
with increasing matter density, and consequently a much stronger magnetic field in 
the stars' core. Nevertheless, the interior magnetic fields of neutron stars are not 



Magnetism in Dense Quark Matter 



3 



directly accessible to observation, so one can only estimate their values with heuris- 
tic methods. Estimates based on macroscopic and microscopic analysis, for nuclear 
ll24l [39l , and quark matter considering both gravitationally bound and self -bound 
stars 1 66 1, have led to maximum fields within the range 10 18 - 10 20 G, depending 
whether the inner medium is formed by neutrons l24l[39l . or quarks [66|. 

For instance, from energy-conservation arguments we can estimate the maximum 
field strength for a quark star. One should expect that the magnetic energy density 
does not exceed the energy density of the self-bound quark matter, which is given 
as the energy density at zero pressure that can have a maximum value equal to that 
of the iron nucleus (roughly 939 MeV). Based on this reasoning, the maximum field 
allowed can be estimated as 

enc enc 

From this result, one notice that the inner field can reach values two orders of mag- 
nitude larger than the estimates done for gravitationally bound stars assuming that 
they contain nuclear matter ll24l [39l ■ 

As we will see in Sec. 8, the magnetic field decreases the inner pressure along 
the field direction of the magnetized system. If the magnetic field that makes such a 
pressure component equal to zero is taken as the maximum value of the inner field 
allowable for a stable gravitational bound star, then a star with a quark matter core 
can have a maximum field - 10 19 - 10 20 G HH, while one with nuclear matter 
can only have fields - 10 18 G E4ll39l . 

Therefore, the investigation of the properties of very dense matter in the presence 
of strong magnetic fields is of interest not just from a fundamental point of view, but 
it could be also closely connected to the physics of strongly magnetized neutron 
stars. 



3 Magnetism in Spin-Zero Color Superconductivity 

An important point to keep in mind in our analysis of the field effects in CS is that 
in spin-zero color superconductivity the electromagnetism is not the conventional 
one. In the color superconducting medium the conventional electromagnetic field is 
not an actual eigenfield, since it is mixed with one of the gluon fields, much like the 
mixing ocurring between the hyper-field and the W-boson in the electroweak model 
in the presence of the Higgs condensate. Thus, even though the original electro- 
magnetic U{\) em symmetry is broken by the formation of the charged quark Cooper 
pairs in the CFL phase J6), a residual t/(l) gauge symmetry still remains. The mass- 
less gauge field associated with this symmetry is given by the linear combination of 
the conventional photon field and the 8''' gluon field Il6ll7l l75l . 

An =cos&4„ -sinQG 8 (2) 
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The corresponding orthogonal linear combination 




sin&4 M + cos0G^ 



(3) 



is massive. The field plays the role of an in-medium or rotated electromagnetic 
field. A magnetic field associated with can penetrate the CS without being sub- 
ject to the Meissner effect, since the color condensate is neutral with respect to the 
rotated charge. However, the rotated electromagnetic field in the CFL superconduc- 
tor is mostly formed by the original photon with only a small admixture of the 8'' ! 
gluon since the mixing angle, cos 6 = g/y / e 2 /3 + g 2 , is sufficiently small. 

The generator of the unbroken U(l) symmetry, which corresponds to the long- 
range rotated photon in the CFL phase, is a matrix in flavor^T,) ® color (3*3) 
space given by Qcfl = Q® I + 1 ® 7g/v3, where Q is the conventional electro- 
magnetic charge operator of quarks and T% is the 8 ,/! Gell-Mann matrix. Using 
the matrix representations, Q — diag(— 1/3, — 1/3,2/3) for (s,d,u) flavors, and 
Tg = diag(— l/\/3, — 1 / a/3, 2/ v^3) for {b,g,r) colors, the Q charges (in units of 
e = ecos 0) of different quarks are 



For the 2SC color superconductor the ground state is formed by spin-zero di- 
quarks, which are also neutral with respect to the rotated electromagnetic field asso- 
ciated with the remnant U(l) symmetry. In this system (see Refs. lB1 [T7l[3Tlll22L 
the generator of the remnant U ( 1 ) symmetry is a matrix in flavor^ x 2) ® color q x 3) 
space given by Qisc = Q®I — I®T% /y/3, with the usual matrix of electromagnetic 
charges of quarks in flavor space Q = diag(2/3, — 1 /3), and 7s is the eighth gener- 
ator of the SU (3) c gauge group in the adjoint representation. The rotated charges of 
the quarks in units of e = ecos 9, are given in this phase by 



and the massless rotated electromagnetic field and orthogonal massive field are de- 
fined in the same way as in the CFL case. 

From now on, we will use "magnetic field" in short, when we refer to the "rotated 
magnetic field," since inside the superconductor only the rotated magnetic field is 
the physical long range field. 



Sb Sg S r dfr dq d r lib u g u r 



(4) 



00-10 0-1+1+10 




(5) 



4 The Magnetic CFL Phase 



The fact that the rotated magnetic field can penetrate the spin-zero color supercon- 
ductor brings the possibility to look for possible field-interaction effects on the CS 
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phase. An important consequence of this interaction was first studied in II561 . It 
is based in the following observation. Although the Cooper pairs have zero rotated 
charge, they can be formed either by neutral quarks or by quarks of opposite charges. 
If the magnetic field is strong enough so that the magnetic length — \[2eB be- 
comes smaller than the pairs' coherence length, then the magnetic field can interact 
with the pair constituents and significantly modify the pair structure of the conden- 
sate. As shown in [56, 57l l58ll59ll . the presence of a magnetic field changes the CFL 
phase characterized by one single gap, producing a splitting of the CFL gap into a 
gap that gets contributions from both pairs of oppositely charged, as well as neutral, 
quarks, denoted by A B , and one that only gets contributions from pairs of neutral 
quarks, denoted by A . The new phase that forms in the presence of the magnetic 
field also has color-flavor-locking, but with a smaller symmetry group SU(2)c+l+r, 
a change that is reflected in the splitting of the A and A B gaps. The less symmetric re- 
alization of the CFL pairing that occurs in the presence of a magnetic field, is known 
as the magnetic-CFL (MCFL) phase |E3|57]|l8]|59]. The MCFL phase has similari- 
ties, but also important differences with the CFL phase l56ll571l58ll59ll62lllllllll2L 
In the strong-field limit, the gap formed by pairs of neutral and charged quarks 
satisfies the gap equation ||56l l57l 1581 l59ll 

j _ g 2 f d 3 q 1 ^ g 2 eB [ A dq 1 



3A 2 J A {In? ^/(q~^ + 2(A B ) 2 6A 2 J- A (2k) 2 ,/ (q ~ n) 2 + W ' 

(6) 

The solution of (|6|l is given by 

/ 3A 2 n 2 \ 

A B ^2y/8fiexp[ , (7) 

V g 2 (^ 2 + f) J 
which can be compared with the CFL gap 

/ 3A 2 n 2 \ 

Acfl ~ 2 exp ^ - j , (8) 

Here we used 8 = A — jj., with A the ultraviolet cutoff of the NJL model that should 
be much larger than any of the typical energy scales of the system, and ji the baryon 
chemical potential. 

The gap A, formed only by pairs of neutral quarks, should be found as the solu- 
tion of the gap equation 

lw jL/- A_flZ 1 + 1 1 ) (9) 

4A2 J A (2%f \ 9 ^(q-pLf-A 2 9 ^(q-pLf + 2(A B ) 2 ) ' 

where it is apparent the interconnection with the gap A B . This is how through A B 
the magnetic field can affect A although it is formed only by neutral quarks as we 
already pointed out. 

The solution of Q is given by 
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1 / 36 21 1 3 \ 

A - exp ( Wx + ffMi+y) + 2-J AcFL ' (10) 

where x = g 2 fx 2 /A 2 7t 2 , andy = eB/jX 2 

The exponent in (|7]i has the typical BCS form, but with^different density of states 
for neutral and charged quarks, i.e. exp [ 1/(A^ +Ng) G ], where = jX 2 /n 2 is 
the density of states at the Fermi surface of the neutral quarks with single chirality, 
Ng = eB/2n 2 is the density of states of the charged quarks lying at the zero Landau 
level at the Fermi surface, and G = —g 2 /3A 2 is the characteristic effective coupling 
constant of the 3 channel 11341 . The effect of the strong magnetic field eB/2 > ji 2 
is to increase the total density of states, thus producing a gap enhancement. 

Although the situation here shares some similarities with the magnetic catalysis 
of chiral symmetry breaking 113|55j|Pl|S]|771|78l|90l|S|^|96l; the way the field 
influences the pairing mechanism in the two cases is quite different. The particles 
participating in the chiral condensate are near the surface of the Dirac sea. The effect 
of a magnetic field there is to effectively reduce the space dimension where the 
particles are embedded at the lowest Landau level (LLL), which as a consequence 
strengthens their effective coupling, and so catalyzing the chiral condensate. Color 
superconductivity, on the other hand, involves quarks near the Fermi surface, with a 
pairing dynamics that is already (1 + 1) -dimensional. Therefore, the B field does not 
yield further dimensional reduction of the pairing dynamics near the Fermi surface 
and hence the LLL does not have a special significance here. Nevertheless, the field 
increases the density of states of the Q-charged quarks, and it is through this effect, 
as shown in Eq. ([7]), that the pairing of the charged particles is reinforced by the 
penetrating magnetic field. 

Note that our analytic solutions are only valid at strong magnetic fields. For fields 
of this order and larger, the Ab gap is larger than Aqfl at the same density values. 
How fast or slow the gaps do it depends very much on the values of the NJL cou- 
plings. For example, for x ~ 0.3, one finds A ~ 0.2 Ag for y = 3/2, while for x ~ 1 
then A ~ 0.5 A B . 

In a recent study 11461 . it was discovered that the MCFL phase actually contains 
one more condensate, which we will call Am- This new condensate is associated with 
the magnetic moment of the Cooper pairs. Physically this is easy to understand. The 
presence of a uniform magnetic field explicitly breaks the spatial rotational sym- 
metry 0(3) to the subgroup 0(2) of rotations about the axis along the field. As 
shown in 11461 . this symmetry reduction has non-trivial consequences for the ground 
state structure of the MCFL superconductor. When one performs the Fierz trans- 
formations in the quark system with both Lorentz and rotational 0(3) symmetries 
explicitly broken, various new pairing channels appear allowing in principle the for- 
mation of new condensates. Of particular interest is an attractive channel that leads 
to a spin-one condensate of Dirac structure Am ~ C/s/'y 2 . Such a gap does not 
break any symmetry that has not already been broken by the other condensates of 
the MCFL ground state, so it in principle is not forbidden. The new condensate cor- 
responds to the zero spin projection of the average magnetic moment of the Cooper 
pairs in the medium. 
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From a physical point of view, it is natural to expect the formation of this extra 
condensate in the magnetized system because the diquarks formed by oppositely 
charged quarks with opposite spins will have a net magnetic moment that may point 
parallel or antiparallel to the magnetic field. Diquarks formed by quarks lying on 
any non-zero Landau level can have magnetic moments pointing in both directions, 
because each quark in the pair may have both spins. Hence the contribution of these 
diquarks to the net magnetic moment should tend to cancel out. On the other hand, 
diquarks from quarks in the LLL can only have one orientation of their magnetic 
moment with respect to the field, because the quarks in the LLL have only one pos- 
sible spin projection. This implies that the main contribution to the new condensate 
should come from the quarks at the LLL, an expectation that is consistent with the 
numerical results found in [46 1, where the new gap was obtained to be negligibly 
small at weak magnetic fields, where the zero Landau level occupation is not sig- 
nificant. On the other hand, at strong magnetic fields, the condensate became com- 
parable in magnitude to the original condensates, A and As, of the MCFL ground 
state [56 1, because the majority of the quarks occupy the LLL in that case. 

Although this new condensate is zero at zero magnetic field, we cannot ignore it 
even at very small magnetic fields because a self-consistent solution of the gap equa- 
tions with A 7^ 0, and Ag ^ 0, but Am = is not possible. This is easy to understand 
since, as long as the magnetic field is not zero, there is always some occupation of 
the LLL. Thus, once a magnetic field is present, Am has to be considered simulta- 
neously with the spin-zero MCFL gaps. 

The Am condensate of the MCFL scenario described above shares a few simi- 
larities with the dynamical generation of an anomalous magnetic moment recently 
found in massless QED l64l l65ll . Akin to the Cooper pairs of oppositely charged 
quarks in the MCFL phase, the fermion and antifermion that pair in massless QED 
also have opposite charges and spins and hence carries a net magnetic moment. A 
dynamical magnetic moment term in the QED Lagrangian does not break any sym- 
metry that has not already been broken by the chiral condensate. Therefore, once the 
chiral condensate is formed due to the magnetic catalysis of chiral symmetry break- 
ing ll54l 15511771 78, 90, 91 96, 97|, the simultaneous formation of a dynamical mass 
and a dynamical magnetic moment is unavoidable 1641 1651 . The realization of the 
anomalous magnetic moment condensate in magnetized massless QED produces a 
non-perturbative Zeeman effect |64[ |65l . 

At moderate magnetic fields the energy gaps A and Ag exhibit oscillations when 
eB/n 2 is varied 1111111121 . owed to the de Haas-van Alphen effect [80, 81] typical 
of charged fermion systems under magnetic fields (see for instance ATI l42l II 111 
112Q, while for Am the oscillations are almost absent [46|. These features indicate, 
as already pointed out, that the main contribution to Am should come from pairs 
whose charged quarks are at the LLL. 

The previous discussion can be visualized in the plot of the gaps as functions of 
a dimensionless parameter eB/ji 2 given in Fig[T] Note that for small magnetic field, 
A and Ag are close to each other and approach the CFL gap Aqfl = 25MeV. As 
the magnetic field increases, A and Ag display oscillatory behaviors with respect 
to eB/jX 2 as long as SB < jx 2 . As originally explained by Landau [93 1, these os- 
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Fig. 1 The three gaps of the MCFL phase as a function of eB/ fl 2 for fl ■ 
with respect to the CFL gap Acfl = 25MeV. 



■ 500MeV. They are scaled 



dilations reveal the quantum nature of the interaction of the charged particles with 
the magnetic field (the well-known Landau quantization phenomenon), and are pro- 
duced by the change in the density of states when passing from one Landau level to 
another. The oscillations cease when the first Landau level exceeds the Fermi sur- 
face. For ultra-strong fields, when only the LLL contributes to the gap equation, Ab 
is much larger than A, as it was found by analytical calculation in l56l . 

The only contribution to Am from higher LLs can come when the number of 
particles is odd, so there are energy states occupied by a single particle, but that is 
a very small part. The cancelation does not occur, however, between the pairs of 
quarks in the LLL because they can only be formed by positive quarks with spin up 
and negative quarks with spin down. At low fields, the number of quarks in the LLL 
is scarce, while for fields of order eB > fl 2 , all the particles are constrained to the 
LLL, hence the variation of Am from lower values at weak field, to higher values at 
sufficiently strong fields. 

It is apparent from the graphical representation of Am in Fig[T[ that its value 
remains relatively small up to magnetic-field values of the order of fx 2 . In the field 
region between 10 18 - 10 19 G, the magnitude of Am grows from a few tenths of Mev 
to tens of Mev. It becomes comparable to the MCFL gap Ab when the field is strong 
enough to put all the quarks in the LLL, shown in the final segment of the plots in 
the figure. 

Another important consequence of the gap Am is the increment in the magnitude 
of Ab for any given value of the magnetic field in the strong field region, as compared 
to its own value found at the same field but ignoring the existence of Am [111, 1 12]. 
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This effect, combined with the increase of Am at strong fields, will make the MCFL 
phase more stable than the regular CFL, a fact that could favor the realization of an 
MCFL core in magnetars. 

Let us now discuss in more detail the difference between the CFL and MCFL 
phases from the point of view of symmetry. In the absence of a magnetic field, 
three-flavor massless quark matter at high baryonic density is in the energetically 
favored CFL phase. There, the diquark condensates lock the color and flavor trans- 
formations, breaking both symmetries. Thus, the symmetry breaking pattern in the 
CFL phase is 

SU(3) C x SU(3) L x SU(3) R x U(l) B -> SU(3) C+L+R . (11) 

In this case, there are only nine Goldstone bosons that survive to the Anderson- 
Higgs mechanism. One is a singlet, scalar mode, associated to the breaking of the 
baryonic symmetry, and the remaining octet is associated to the axial SU(3)a group, 
just like the octet of mesons in vacuum. At sufficiently high density, the anomaly 
is suppressed, and then one can as well consider the spontaneous breaking of an 
approximated U{\)a symmetry, and the additional pseudo Goldstone boson. We 
will ignore this effect, though. 

Once electromagnetic effects are considered, the flavor symmetries of QCD are 
reduced, as only the d and s quarks have equal electromagnetic charges, q = —e/3, 
while the u quark has electromagnetic charge, q = 2e/3. However, because the elec- 
tromagnetic structure constant a e . m . is relatively small, this effect is considered to 
be really tiny, a small perturbation, and one can consider good approximated fla- 
vor symmetries. Nevertheless, in the presence of a strong magnetic field one cannot 
consider the effects of electromagnetism as a small perturbation. Flavor symmetries 
are explicitly reduced from SU(3)l,r to SU(2) LR . For sufficiently strong magnetic 
fields, in the MCFL phase the symmetry breaking pattern is then 

SU(3) C x SU(2) L x SU(2) R X U(l) A l) x U(l) B x C/(l) e . m . -> SU{2) c+l +r- (12) 

Here the symmetry group U(l) is related to a current, which is an anomaly free 

linear combination of u,d and s axial currents, and such that U(l)^ C SU(3)a- The 
locked SU (2) group corresponds to the maximal unbroken symmetry, such that it 
maximizes the condensation energy. The counting of broken generators, after taking 
into account the Anderson-Higgs mechanism, leads to only five Goldstone bosons. 
As in the CFL case, one is associated with the breaking of the baryon symmetry; 
three Goldstone bosons are associated with the breaking of SU (2) a, and another one 
with the breaking of U (l)j^. As before, if the effects of the anomaly could be ne- 
glected, there would be another pseudo Goldstone boson associated with the U(1)a 
symmetry. Thus, apart from modifying the value of the gaps, an applied strong mag- 
netic field also affects the number of Goldstone bosons, reducing them from nine 
(neutral and charged) to five (neutral). 

Once a magnetic field is present, the original symmetry group is reduced, and the 



low energy theory correspond to the breaking pattern ( 12 1, hence be described by 
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five Goldstone bosons. In practice however, at weak magnetic fields, it is reasonable 
to treat the symmetry of the CFL phase as a good approximated symmetry, which 
means that at weak fields the low-energy excitations are essentially governed by 
nine approximately massless scalars (those of the breaking pattern (Hi) instead of 
five. 

A question of order here is: what do we exactly understand as a weak magnetic 
field? In other words, what is the threshold-field strength that effectively separates 
the CFL low energy behavior from the MCFL one? A fundamental clue in this direc- 
tion was found in l62l by determining the term in the low-energy CFL Lagrangian 
that generates a field-induced mass for the charged Goldstone fields, so disconnect- 
ing them from the low-energy dynamics at some field strength and thereby effec- 
tively reducing the number of Goldstone bosons from the nine of the CFL phase, to 
the five neutral ones of the MCFL. 

The threshold field Bmcfl for the effective CFL — > MCFL symmetry crossover 
was found to be [62 1 

eBMCFL = ^^CFL - ^CFL- (13) 

where the weak-field approximation vj_ ~ 1 /y3 was considered H135L The thresh- 
old field does not depend on the decay constant f x , therefore it depends on ll 
only through Acfl- For Acfl ~ l5MeV one gets SBmcfl ~ 10 16 G. At these field 
strengths, the charged mesons decouple from the low-energy theory. When this de- 
coupling occurs, the five neutral Goldstone bosons (including the one associated 
to the baryon symmetry breaking) that characterize the MCFL phase will drive the 
low-energy physics of the system. Therefore, coming from low to higher fields, 
the first magnetic phase that will effectively show up in the magnetized system for 
eB ~ A CFL will be the MCFL ED EZ1 E3 ED . 

Summarizing, in a color superconductor with three-flavor quarks at very high 
densities an increasing magnetic field produces a phase crossover from CFL to 
MCFL. During this phase transmutation no symmetry breaking occurs, since in 
principle once a magnetic field is present the symmetry is strictly speaking that 
of the MCFL, as discussed above. However, in practice for B ~ Bmcfl ~ A CFL trie 
main features of MCFL emerge through the low-energy behavior of the system l62l . 
At the threshold field Bmcfl, only five neutral Goldstone bosons remain out of the 
original nine characterizing the low-energy behavior of the CFL phase, because the 
charged Goldstone bosons acquire field dependent masses and can decay in lighter 
modes. For a meson to be stable in this system, its mass should be less than twice 
the gap, otherwise it will decay into a particle-antiparticle pair. That means that, as 
proved in Ref. [62], once the applied field produces a mass for the charged Gold- 
stones of the order of the CFL gap it is reached the threshold field for the effective 
CFL — > MCFL symmetry transmutation. 

The existence of this phase transmutation is on the other hand manifested in the 
behavior of the gaps versus the magnetic field. At field strength smaller than the 
threshold field we find that A s» Ag ~ Acfl, while for fields closer to Bmcfl the 
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gaps exhibit oscillations with respect to SB/ jx 1 Il46lllllllll21 . owed to the de Haas- 
van Alphen effect ||80ll8fl . 

It is worth to call attention to the analogy between the CFL-MCFL crossover and 
what could be called a field-induced Mott transition. Mott transitions were origi- 
nally considered in condensed matter in the context of metal-insulator transitions 
in strongly-correlated systems II 1071 1 1081 . Later on, Mott transitions have been also 
discussed in QCD to describe derealization of bound states into their constituents 
at a temperature defined as the Mott temperature [82 1. By definition, the Mott tem- 
perature Tm is the temperature at which the mass of the bound state equals the mass 
of its constituents, so the bound state becomes a resonance at T > Tm- In the present 
work, the role of the Mott temperature is played by the threshold field Bmcfl- Mott 
transitions typically lead to the appearance of singularities at T = Tm in a num- 
ber of physically relevant observables. It is an open question, worth to be investi- 
gated whether similar singularities are or not present in the CFL-MCFL crossover 
at Bmcfl- 



5 Magnetoelectric Effect in Cold-Dense Matter 

It is well known that the phenomenon of CS shares many characteristics of con- 
densed matter systems lfT2l . In this section, we discuss a new feature of CS that has 
its counterpart in magnetically ordered materials and has been known in the context 
of condensed matter for many years. It is the so called magnetoelectric (ME) effect, 
which establishes a relation between the electric and magnetic properties of certain 
materials. In general, it states that the electric polarization of such materials may de- 
pend on an applied magnetic field and/or that the magnetization may depend on an 
applied electric field. The first observations of magnetoelectricity took place when 
a moving dielectric was found to become polarized when placed in a magnetic field 
11125111411 . In 1894, Pierre Curie 1 38 1 was the first in pointing out the possibility of 
an intrinsic ME effect for certain (non-moving) crystals on the basis of symmetry 
considerations. But it took many decades to be understood and proposed by Lan- 
dau and Lifshitz [95| that the linear ME effect is only allowed in time-asymmetric 
systems. Recently the ME effect regained new interest in condensed matter thanks 
to new advancements in material science and with the development of the so-called 
multiferroic materials for which the ME effect is significant for practical applica- 
tions ED. 

As demonstrated in Refs. [47, 48] the ME effect also occurs in a highly magne- 
tized CS medium like the MCFL phase. In particular, in ll47ll48ll it was shown how 
the electric susceptibility of this medium depends on an applied strong magnetic 
field. 

Let us start by discussing the ME effect at weak fields. At weak fields this effect 
can be studied by taking into account the expansion of the system's free energy in 
powers of the electric E and magnetic B fields 
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F(E, B) = F - OdEi - PiBt - jijEiBj - ^jEiEj - z^Bfij 

— KijkEiEjB k — XijkEiB jB k — OjjkiEjE jB k B[ — ... ( 14) 

In this weak-field expansion the coefficients Of;, ju etc., which are the suscepti- 
bility tensors, can be found from the infinite set of one-loop polarization operator 
diagrams with external legs of the in-medium photon field and internal lines of 
the full CFL quark propagator of the rotated charged quarks. Hence, these coeffi- 
cients can only depend on the baryonic chemical potential, the temperature and the 



CFL gap. From ( 14 1, the electric polarization can be found as 



dF ~ ~ ~~ „ ~ ~ 

Pi = --~= a, + y ij B j + 2ri ij E j + 2K i j k EjB k + X ijk B j B k + 2a ijk iE j B k B l + .. . (15) 
dEi 

If the tensor y is different from zero the system exhibits the linear ME effect. 



From the free energy ( 14 1, we see that the linear ME effect can only exist if the 
time-reversal and parity symmetries are broken in the medium. In the CFL phase, 
the time-reversal symmetry is broken by the CFL gap [46 1, but parity is preserved. 
Thus, the linear ME effect cannot be present in this medium. The behavior under a 
time-reversal transformation underscores an important difference between the CFL 
color superconductivity and the conventional, electric superconductivity. While the 
CFL color superconductor is not invariant under time-reversal symmetry, the con- 
ventional superconductor is, since in the conventional superconductor the Cooper 
pairs are usually formed by time-reversed one-particle states [16]. In the conven- 
tional superconductor the violation of the T-invariance occurs only via some exter- 
nal perturbation which can lead in turn to pair breaking and to the so-called gapless 
superconductivity (2). 

Higher-order ME terms are parameterized by the tensors fc, and X . As it happens 
with 7, the coefficient X ^ is forbidden because it requires parity violation. On 
the other hand, although a K ^ term only requires time-reversal violation, to form 
a third-rank tensor independent of the momentum and parity invariant, the medium 
would need to have an extra spatial vector structure. However, the only tensor struc- 
tures available to form such a third-rank tensor in the CFL phase are the metric 
tensor g^ v and the medium fourth velocity , which in the rest frame is a temporal 
vector 11^ = (1,0,0,0), so the coefficient K should be zero too. Hence, we do not 
expect any ME effect associated with the lower terms in the weak-field expansion 



of the free energy ( 15 i 



At strong magnetic fields, the situation is quite different. In this case the expan- 
sion of the free energy can only be done in powers of a weak electric field, and the 
coefficients of each term can be found from the corresponding one-loop polarization 
operators, which now depend on the strong magnetic field in the MCFL phase. The 
free energy expansion in this case takes the form 

F' (E,B) = (5) - afc - TjljEiEj (16) 
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The tensors a' and T]' can depend now on the baryonic chemical potential, temper- 
ature, magnetic field and gaps of the MCFL phase [46|. They can be found respec- 
tively by calculating the tadpole and the second rank polarization operator tensor of 
the MCFL phase in the strong field limit. An a' ^ would indicate that the MCFL 
medium behaves as a ferroelectric material [86 , 99], but this is not the case because 
this phase is parity symmetric R6ll . hence a' = 0. The tensor r\' ', nevertheless, is 
not forbidden by any symmetry argument. If it is different from zero, r\' would 
characterize the lowest order of the system dielectric response. More important, if 
r\' results to be dependent on the magnetic field, this would imply that the electric 
polarization P = f]'E depends on the magnetic field through r\' , hence the MCFL 
phase would exhibit the ME effect. 

To find the electric susceptibility 77/ in the strong-magnetic -field limit of the 
MCFL phase we start from ll46l 

F'(E,B) -^(B) ~ - J A (x3)n 00 {x 3 ~ x' 3 )A (x' 3 )dx 3 dx 3 = -ri'E 2 , (17) 

our task is then reduced to the calculation of the zero-zero component of the one- 
loop polarization operator at strong magnetic field in the infrared limit, Iloo(po = 
Q,p->0). 

Now, the photon polarization operator should be gauge invariant. That is, in the 
strong-field approximation, it should satisfy the transversality condition in the re- 
duced (1 + 1)-D space (Pn"ll]iv(p" ) = 0)- As known, the polarization operator ten- 
sor can be expanded in a superposition of independent transverse Lorentz tensors. 
The number of these basic transverse tensors depends on the symmetries of the sys- 
tem under consideration. For example, in vacuum, where the only available tensorial 
structures are the four-momentum and the metric tensor, there is only one gauge in- 
variant structure. When a medium is under consideration (i.e. at finite temperature 
or finite density), since the Lorentz symmetry is broken, there is an additional gauge 
invariant structure that can be formed by taking into account a new four-vector, the 
four- velocity of the medium center of mass, u^, [68 1. When a magnetic field is ap- 
plied on that medium, then the structure of the polarization operator is enriched 
by an additional tensor, F^ v . Then, at finite density and in the presence of a mag- 
netic field, there are nine independent gauge-invariant tensorial structures 11 161 . At 
strong magnetic field, when the particles are confined to the LLL, due to the fact 
that the transverse momentum is zero, there is a dimensional reduction leaving only 

the tensors g^ v ,p^ and a| = (1,0) at our disposal. The original nine structures in 
[ 1 16 1 now reduce to only two 

(1) 11 9 II 

7/Iv = (PV gfiv-pjiPv, ( 18 ) 

and 

1_pU^1], (19) 

p 1 



T [Z) _ 



(P 1 



■P" 

\2 
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Moreover, one can readily check that the two tensors ( |T8] > and ( [T9| > are equivalent, 
which indicates that the rotated-photon polarization operator tensor, at strong mag- 
netic field, only has one independent structure 



njivip 11 ) = n(p\n,B) (/) g^v-pfiPv 



(20) 



with n(p>> ,n,B) being a scalar coefficient depending on the photon longitudinal 
momentum, baryonic chemical potential and magnetic field. 

At zero temperature, the regularized components of the polarization operator in 
powers of the photon momentum components po and p$, up to quadratic terms are 
given by 



^2 1^/5 



- lim 

A-t«> 



eB\pj 
67t 2 



1 

a 2 



- lim 

A- 



e 2 \eB\ 
6n 2 



2 

Po 
A 2 



1 

A 2 



eB\ P j 



6k 2 A 2 



2 

Po 
A 2 



eB\p 2 



6k 2 A 2 



(21) 



(22) 



and JJj,or = TIqt, r ~ 0. As should be expected, the regulator A introduced through 
the Pauli- Villars regularization scheme does not appear in the final results once we 
take A^oo. 

Because Uoo has no constant contribution in the infrared limit po = 0,p3 —> 0, one 
immediately concludes that there is no Debye screening in the strong-field region, 
as it was the case at zero field in the CFL phase 1 100 , 127|. This is simply because 
all quarks are bound within the rotated-charge neutral condensates. There is also no 
Meissner screening (i.e. is zero in the zero-momentum limit), as it should be 
expected from the remnant U(l) gauge symmetry. However, the condensates have 
electric dipole moments and could align themselves in an electric field. Hence, this 
should modify the dielectric constant of the medium. Since the quadratic term in the 
effective U (1) Lagrangian is given by A^(— p)[D~l(p) +n flv (p)]A v (p), with D 1 
being the bare rotated photon propagator, the effective action of the 0(1) field in the 
strong-field region is given by 



2 ± 



E 



2An 11 



H 



1 

21, 



H ± H 



(23) 



where the separation between transverse and longitudinal parts is due to the 0(3) —> 
0(2) symmetry breaking produced by the strong magnetic field B. In (23 i, E, H 



are weak electric and magnetic field probes, respectively. In ( 23 1 the coefficients 



e and A denote the electric permittivity and magnetic permeability of the medium 
respectively. 



From (21 



-( 22 1 it is straightforward that in the infrared limit the transverse and 
longitudinal components of the electric permittivity and magnetic permeability be- 
come 

* 2 \eB\ 



A ± =A||~1, £ ± = 1, £|| = 1+Xmcfl = 1 



6k 2 A 2 ' 



(24) 
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where Xmcfl * s tne longitudinal electric susceptibility. Notice that the longitudinal 
electric susceptibility is much larger than one because in the strong-magnetic-field 
limit eB->Al ll46l . 

Although a static U(l) charge cannot be completely Debye screened by the 0(1) 
neutral Cooper pairs, it can still be partially screened along the magnetic field di- 
rection because the medium is highly polarizable on that direction. This is due to 
the existence of Cooper pairs with opposite rotated charges Q that behave as electric 
dipoles with respect to the rotated electromagnetism of the MCFL phase. Moreover, 
the electric susceptibility depends on the magnetic field. When the magnetic field 
increases in the strong-field region, the susceptibility becomes smaller, because the 
coherence length % ~ l/4o decreases (i.e. Aq increases) with the field at a quicker 
rate than VsB [46], and the pair's coherence length £, plays the role of the dipole 
length. Hence, with increasing magnetic field the polarization effects weaken in the 
strong-field region. The tuning of the electric polarization by a magnetic field is what 
is called in condensed matter physics the magnetoelectric effect. From ( 24 1, we also 
see that at strong magnetic fields the medium turns out to be very anisotropic. The 
fact that the electric permittivity is only modified in the longitudinal direction is due 
to the confinement of the quarks to the LLL at high enough fields. 



6 Paramagnetism in Color Superconductivity 

Another nontrivial electromagnetic effect in cold-dense QCD is that an applied mag- 
netic field can interact inside the color superconductor with the gluons, which as 
known, are neutral with respect to the conventional electromagnetism in vacuum. 

Thus, we now analyze how the gluons are affected by an applied magnetic field in 
a CS state and how at sufficiently strong magnetic fields a new phase, that we call the 
Paramagnetic-CFL (PCFL) phase 116011611 . is created. In the color superconductor 
some of the gluons acquire rotated electric charges. In the CFL phase the g-charge 
of the gluons in units of e are 



G n G n G fi G ti G fi l \i 1 \i G n 
1 -1 1 -1 



(25) 



The Q-charged fields in (25 1 correspond to the combinations = A= [G 4 ^ =F iG 5 ^} 

m&I± = ±[GlTiGl\. 

To investigate the effect of the applied rotated magnetic field H on the charged 
gluons, we should start from the effective action of the charged fields G^ (the con- 
tribution of the other charges gluons is similar) 

r eff = J dx{-^(f^ v ) 2 + G M [(n Al n |i )5 MV -2ief^ v 
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- (modpoSvo+mlfSpidvi) - (1 - - J7 M JT V )]G+} (26) 

Here, g is the gauge fixing parameter, = <9^ — ieA^ is the covariant derivative in 
the presence of the external rotated field, mo and rtiM are the G^ -field Debye and 
Meissner masses respectively, and the field strength tensor for the rotated electro- 
magnetic field if denoted by fay — d^Ay — d v A^. The corresponding Debye and 
Meissner masses in (26 1 are given by 1100111271 



2 2 21 -81n2 2 2 21 -81n2 
D = g 18 ' ™ M = 8 54 ' 

with m\ =g 2 {jx 2 /2n 2 ). We are neglecting the correction produced by the applied 
field to the gluon Meissner masses since it will be a second order effect. The ef- 



fective action (26 1 is characteristic of a spin-1 charged field in a magnetic field (for 
details see for instance P31I441 ). 

Assuming an applied magnetic field along the third spatial direction { fff = H), 
we find after diagonalizing the mass matrix of the field components (G\ , G^) in 
& 

m^ieH\fml+m \ 
-ieH nty J \ mjj — eH J 

with corresponding eigenvectors (Gj ,Gi~) — ► (G,iG). We see that the lowest mass 



mode in ( 28 1 has a sort of "Higgs mass" above the critical field eHc = »%, indicat- 
ing the setup of an instability for the G-field. This phenomenon is the well known 
"zero-mode problem" found in the presence of a magnetic field for Yang-Mills fields 
II 1 101 1 1321 . for the bosons in the electroweak theory M [TBI [1331 . and even for 



higher-spin fields in the context of string theories |49 53 1 and it is due to the pres 



ence of the gluon anomalous magnetic moment term 2/e/^ v G^Gy in (26 1. Thus, to 
remove the instability it is needed the restructuring of the ground statethrough the 
condensate of the field bearing the tachyonic mode (i.e. the G-field). 

It is worth to call attention that the gluon condensate under consideration is not 
the only charged spin-one condensate generated in a theory with a large fermion 
density. As known B0ll5Tll98l . a spin-one condensate of W-bosons can be origi- 
nated at sufficiently high fermion density in the context of the electroweak theory 
at zero magnetic field. However, the physical implications of the gluon condensate 
induced by the magnetic field in the CS are fundamentally different from those as- 
sociated to the homogeneous W-boson condensate of the dense electroweak theory 
ll50l ISTl 198 1 . One of the main physical differences is that the homogeneous W con- 
densate, being electrically charged, so to compensate the excess of charge due to the 
finite density of electrons [50, 51 , 98 1, breaks the electromagnetic U(l) group pro- 
ducing a conventional superconducting state ||52l ; while the inhomogeneous gluon 
condensate in CS is formed with gluons of both charges, so keeping the condensate 
state neutral. 

To find the G-field condensate and the induced magnetic field B = V x A, with 
A being the total rotated electromagnetic potential in the condensed phase in the 
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presence of the external field H, we should start from the Gibbs free energy density 
Sf = & — HB, since it depends on both B and H (JF is the system free energy 
density). Since specializing H in the third direction the instability develops in the 
(x,y) -plane, we make the ansatz for the condensed field G — G(x,y). Starting from 



(26 1 in the Feynman gauge g = 1, which in terms of the condensed field G implies 



(iTi + iTl2)G = 0, we have that the Gibbs free energy in the condensed phase is 

% = &n0 + n 2 G 2 -2{eB- m 2 M )G 2 + 2g 2 G 4 + l -B 2 - HB. (29) 

where J£"„o is the system free energy in the normal phase (G = 0) at zero magnetic 
field. 

The minimum equations for the fields G and B are respectively obtained from 
([29) as 

n 2 G + 2{m 2 M -eB)G + 8g 2 G 2 G = 0, (30) 
2eG 2 -B + H = (31) 



Identifying G with the complex order parameter, Eqs. (30i-(31 1 become analogous 



to the Ginzburg-Landau equations for a conventional superconductor except by the 



negative sign in front of the B field in Eq. ( 30 1 and the positive sign in the first term 



of the LHS of Eq. (31 1 |60|. The fact that those signs turn the opposite of those 
appearing in conventional superconductivity is due to the different nature of the 
condensates in both cases. While in conventional superconductivity the Cooper pair 
is a spin-zero condensate, here we have a condensate formed by spin-one charged 
particles interacting through their anomalous magnetic moment with the magnetic 



field (i.e. the term 2ief flv G ll Gy in (26 1) 



Notice that because of the different sign in the first term of ( 3 1 1, the resultant 
field B is stronger than the applied field H, contrary to what occurs in conventional 
superconductivity. Thus, when a gluon condensate develops, the magnetic field will 
be antiscreened and the color superconductor will behave as a paramagnet. The an- 
tiscreening of a magnetic field has been also found in the context of the electroweak 
theory for magnetic fields H > Mm/e ~ 10 24 G lfT4l . Just as in the electroweak case, 
the antiscreening in the color superconductor is a direct consequence of the asymp- 
totic freedom of the underlying theory lfl4l . 

Therefore, the magnetic field in the new phase is boosted to a higher value, which 
depends on the modulus of the G-condensate. That is why the phase attained at 
H > H c is called paramagnetic CFL (PCFL) EHUD. It should be pointed out that 
at the scale of baryon densities typical of neutron-star cores (ji ~ AQQMeV , g(ll) — 
3) the charged gluons magnetic mass in the CFL phase is ~ 16 x 10~ 3 GeV 2 . 
This implies a critical magnetic field of order H c ~ 0.7 x 10 17 G. Although it is 
a significant high value, it is in the expected range for the neutron star interiors 
with cold-dense quark matter ||66l 11131 . Let us underline that in our analysis we 
considered asymptotic densities where quark masses can be neglected. At lower 
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densities where the Meissner masses of the charged gluons become smaller, the 
field values needed to develop the magnetic instability will be smaller. 

To find the structure of the gluon condensate we should solve the non-linear dif- 



ferential equation (30i. However, to get an analytic solution we can consider the 
approximation where H m H c — m 2 ^ and consequently | G | « 0. In this approxima- 
tion, Eq. ( [30] ) can be linearized as 

rV3 A%i~ - , B 2 , 1 ._ 

[df - ^Bxdy - An 2 -^-x 2 -—]G = 0, j=x,y (32) 
O <2>o 

where we fixed the gauge condition A 2 — Bx\, and introduced the notations <J>o = 
2n/e, and£ 2 = \(eB ~m 2 M )- 1 . 

Eq. < [32] > is formally similar to the Abrikosov's equation in type-II conventional 
superconductivity El 151 1751 , with % playing the role of the coherence length and <Pq 
of the flux quantum per vortex cell. Then, following the Abrikosov's approach, a 



solution of Eq. ( 32 1 can be found as 



G{x. /y ) = -j=—e ^(ii/t), (33) 



with d-i(u/z) being the elliptic theta function with arguments 



- iKb ^ 2+ h l x= - in Y 2 (34) 



In ( 34 1 the parameter b is the periodic length in the y-direction (b = Ay). The double 
periodicity of the elliptic theta function also implies that there is a periodicity in the 
x-direction given by Ax = <Po/bH c . Therefore, the magnetic flux through each peri- 
odicity cell (AxAy) in the vortex lattice is quantized H c AxAy = <Pq, with <Pq being 
the flux quantum per unit vortex cell. In this semi-qualitative analysis we considered 
the Abrikosov's ansatz of a rectangular lattice, but the lattice configuration should 
be carefully determined from a minimal energy analysis. For the rectangular lattice, 
we see that the area of the unit cell is A — AxAy — &o/H c , so decreasing with H. 

In conclusion, to remove the instability created by an external uniform magnetic 
field in the z-direction, a periodic arrangement of vortices of charged gluon con- 
densates is generated in the (x,y)-plane. The currents in the (^,y)-plane created by 
these vortices increase the magnitude of the net magnetic field in the direction of the 
original field, but since the magnitude of the resultant field varies in the (x,y) -plane, 
the vortex condensate leads to a net inhomogeneous magnetic field. Therefore, the 
presence of a supercritical magnetic field leads to the formation of a fluxoid along 
the z-direction and the appearance of a nontrivial topology on the perpendicular 



plane. From (31 1 we see that the resultant magnetic field can go from a minimum 
value H to a maximum at the core of the fluxoid that depends on the amplitude of 
the gluon condensate determined by the mismatch between the applied field and the 
gluon Meissner mass. 
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Summarizing, at low H field, the CFL phase behaves as an insulator, and the H 
field just penetrates through it without any change of strength. At sufficiently high 
field eH ~ m 2 M , the condensation of G* is triggered inducing the formation of a 
lattice of magnetic flux tubes that breaks the translational and remaining rotational 
symmetries, creating the so called paramagnetic phase. We stress that contrary to 
the situation in conventional type-II superconductivity, where the applied field only 
penetrates through the flux tubes and with a smaller strength, the vortex state in the 
color superconductor has the peculiarity that outside the flux tube the applied field 
H totally penetrates the sample, while inside the tubes the magnetic field becomes 
larger than H (this is the origin of the paramagnetic behavior of this CS phase). This 
effect provides an internal mechanism to increase the magnetic field of a compact 
star with a CS core. 



7 Magnetic Phases in CFL Matter 

From the discussions in the previous sections it is clear that in the three-flavor 
color superconductor at very high densities an increasing magnetic field produces 
a crossover from CFL to MCFL first, and then a phase transition from MCFL to 
PCFL. 

During the crossover, no symmetry breaking occurs, since in principle once a 
magnetic field is present the symmetry is already that of the MCFL. At very weak 
magnetic fields, the color superconducting state is practically described by the CFL 
phase, because the charged mesons corresponding to the Goldstone modes, although 
massive, are so light that they cannot decay in pairs of quark-antiquark. When the 
field strength is of the order of the quarks' energy gap Acfl, the charged mesons 
become heavy enough to decouple and the low-energy physics is indeed that of the 
MCFL phase, where five neutral massless mesons drive the low-energy behavior. 

Going from MCFL to PCFL is, on the other hand, a real phase transition 1601I6TI . 
as the translational symmetry, as well as the remaining rotational symmetry in the 
plane perpendicular to the applied magnetic field are broken by the vortex state. This 
phase transition is driven by fields whose strengths are comparable to the magnetic 
masses mjj of the charged gluons, so creating a chromomagnetic instability that 
leads to the formation of a vortex state and the antiscreening of the magnetic field 
MM- 

This magnetic instability is characteristic of systems of charged bosons with 
higher spins (s > 1). Taking into account that at zero momentum the energy spec- 
trum in a magnetic field H of a charged boson of spin s, charge e, gyromagnetic ratio 
g, and mass m is 

E 2 = (2n + l)eH-geH-s + m 2 , (35) 

it is evident that for spin-one particles, for which g = 2, the energy becomes imag- 
inary, i.e. E 2 < 0, if the field satisfies H > H cr — m 2 /e), implying that when the 
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field surpasses the critical value H cr , one of the modes of the charged gauge field 
becomes tachyonic inducing the vortex formation. 

Within a NJL model, for fields comparable to the baryon chemical potential, the 
ground state is that of the MCFL phase with sizable values of the three condensates 
Ab, Am, and A. However, once the gluon effects are taken into account, the PCFL 
vortex state generated at lower fields is unavoidable and the picture becomes much 
more complicated due to the inhomogeneities of the gluon condensate and net mag- 
netic field. From a physical point of view, it is natural to expect that in this situation 
the three fermion gaps will remain, because their physical origin, is still the same. 
That is, an inhomogeneous magnetic field will also distinguish between pairs of op- 
posite charged quarks and pairs of neutral quarks, and those of opposite charged 
quarks will still have a magnetic moment contributing to the condensate Am- How- 
ever, all these condensates should become inhomogeneous in the (^,y)-plane. 



8 Equation of State of the MCFL Phase 

At present, some of the best-known characteristics of stellar objects are their masses 
and radii. The relation between the mass and the radius of a star is determined by 
the equation of state (EoS) of the inner phase of the matter in the star. If one can 
identify some features connecting the star's internal state (nuclear, strange, color 
superconducting, etc.) to its mass/radius relation, one would have an observational 
tool to discriminate among the actual realization of different star inner phases in 
nature. From previous theoretical studies \M ED US EH HDJ \WT\ [[02] [TT7l 
118 1 26 1 29 ] the mass-radius relationship predicted for neutron stars with different 
quark-matter phases (CS or unpaired) at the core are very similar to those having 
hadronic phases, at least for the observed mass/radius range. As a consequence, it is 
very difficult to find a clear observational signature that can distinguish among them. 
Nevertheless, an important ingredient was ignored in these studies: the magnetic 
field, which in some compact stars could reach very high values in the inner regions. 

As pointed out in [66], a strong magnetic field can create a significant anisotropy 
in the longitudinal and transverse pressures. One would expect then, that the EoS, 
and consequently, the mass-radius ratio, become affected by sufficiently strong core 
fields. Given that we are beginning to obtain real observational constraints on the 
EoS of neutron stars lfl3l . it is important to investigate the EoS in the presence of a 
magnetic field for different inner star phases to be able to discard those that do not 
agree with observations. 

In order to understand the relevance of the magnetic field to tell apart neutron 
stars from stars with paired quark matter, it is convenient to recall that when the 
pressure exerted by the central matter density of neutron stars (which is about 
200 — 600M<?V //m 3 ) is contrasted with that exerted by an electromagnetic field, the 
field strength needed for these two contributions to be of comparable order results 
of order ~ 10 18 G (25). It is worth to notice that even these very strong fields are not 
enough to produce quantum effects like the Landau quantization of the protons, be- 
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cause these effects only show up when the particles's cyclotron energy ehB/mc be- 
comes comparable to its rest energy mc 2 , which for protons means a field ~ 10 20 G. 

However, for stars with paired quark matter, the situation is rather different. 
Naively, one might think that comparable matter and field pressures in this case 
would occur only at much larger fields, since the quark matter can only exist at 
even larger densities to ensure deconfinement. In reality, though, the situation is 
more subtle. As argued in [8 1, the leading term in the matter pressure coming from 
the contribution of the particles in the Fermi sea, ~ fi 4 , could be (almost) canceled 
out by the negative pressure of the bag constant and in such a case, the next-to- 
leading term would play a more relevant role than initially expected. Consequently, 
the magnetic pressure might only need to be of the order of that produced by the par- 
ticles close to the Fermi surface, which becomes the next-to-leading contribution, 
~ H 2 A 2 , with A the superconducting gap and jx the baryonic chemical potential. For 
typical values of these parameters in paired quark matter one obtains a field strength 
~ 10 18 G. Moreover, the magnetic field can affect the pressure in a less obvious way 
too, since as shown in I57ll58ll59l . it modifies the structure and magnitude of the 
superconductor's gap, an effect that, as found in [60, 61 1, starts to become relevant 
already at fields of order 10 16 G and leads to de Haas van-Alphen oscillations of the 
gap magnitude 1111111121 . It is therefore quite plausible that the effects of moder- 
ately strong magnetic fields in the EoS of compact stars with color superconducting 
matter will be more noticeable than in stars made up only of nucleons, where quan- 
tum effects starts to be significant for field four orders of magnitude larger. This is 
why an evaluation of the EoS in magnetized quark phases is necessary and relevant. 

In II 131 . a self-consistent analysis of the EoS of MCFL matter, was performed 
taking into consideration the solution of the gap equations and the anisotropy of the 
pressures in a magnetic field. In that study a uniform and constant magnetic field was 
assumed. The reliability of this assumption for neutron stars, where the magnetic 
field strength is expected to vary from the core to the surface in several orders, is 
based on the fact that the scale of the field variation in the stellar medium is much 
larger than the microscopic magnetic scale for both weak and strong magnetic fields 
l25l . Hence, when investigating the field effects in the EoS, it is consistent to take a 
magnetic field that is locally constant and uniform. This is the reason why such an 
approximation has been systematically used in all the previous works on magnetized 
nuclear Q]|20l|22|28l|2g[3!|35]|7^ and quark 

matter rf32l 1331 l74l H04l H05l H20l . 



8.1 Covariant Structure of the Energy -Momentum Tensor in a 
Magnetized System 

In the reference frame comoving with the many-particle system, the system normal 
stresses (pressures) can be obtained from the diagonal spatial components of the 
average energy-momentum tensor (t' ! ); the system energy, from its zeroth diagonal 
component (t 00 ); and the shear stresses (which are absent for the case of a uniform 
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magnetic field) from the off-diagonal spatial components (t ij ) ll94l . Then, to find the 
energy density and pressures of the dense magnetized system we need to calculate 
the quantum-statistical averages of the corresponding components of the energy- 
momentum tensor of the fermion system in the presence of a magnetic field. 

These calculations were carried out long time ago in Ref. E71 . using a QFT 
second-quantization approach. There, a quantum-mechanical average of the energy- 
momentum tensor in the eigen-states of the Dirac equation in the presence of the 
uniform magnetic field was first performed to get the corresponding quantum op- 
erator in the occupation-number space. The macroscopic stress-energy tensor was 
then found by averaging its quantum operator in the statistical ensemble using the 
many-particle density matrix. Similar calculations were performed in Ref. |66|, but 
using a functional-method approach that makes it easier to recognize the thermo- 
dynamical quantities entering in the final results. An advantage of the procedure 
followed in [66 1, as compared with that of l27l . is that it does not assume that the 
fermion fields entering in the definitions of the energy and pressures satisfy the clas- 
sical equation of motions (i.e. the Dirac equations for y/ and 1/7), but the functional 
integrals integrate in all field configurations. Hence, the terms depending on the La- 
grangian density Jzf v in x^y were kept, while in Ref. |27 1 the condition 5£y= was 
considered as a constraint. 

Let us then consider a hot and dense system of fermions in a constant and uniform 
magnetic field B. At this point it is convenient to introduce the covariant decompo- 
sition for the energy momentum tensor of the whole system containing the matter 
and field contributions. In order to accomplish this goal, we define the system ther- 
modynamic potential as the sum of the matter, £2f, and field, B 2 /2, contributions 

B 2 

n=£2 f +— (36) 

Taking into account the symmetries of the magnetized dense system, we can 
write the statistical average of the energy-momentum tensor as a combination of all 
the available independent structures 

-^-(^ v ) =Qri^ v + {ixN+TS)u^u v +BMri^\ (37) 

where N = — (d£2 /dji)is the particle number density, S — —(d£2/dT^ is the system 
entropy, M = —(dQ/dB) is the system magnetization and J]^ = F^ p F p v (where 
jtW = pVP /g denotes the normalized electromagnetic strength tensor). 

To understand the origin of the covariant decomposition ( |3"7) , notice that as a con- 
sequence of the breaking of the rotational symmetry 0(3) produced by the external 
magnetic field, the Minkowskian metric splits in transverse rj^ v and longitudinal 
T]^ = ri^ v — F^ p Fp structures. Considering the quantum field limit with no mag- 
netic field, i.e. when T = ji = B = 0, the only term different from zero is the first one 
in the RHS of | |3"7) . In that case the system has Lorentz symmetry and the energy 
density, e, and pressure, p, are given by e = — p = £2f. If temperature and/or den- 
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sity are switched on, then the Lorentz symmetry is broken specializing a particular 
reference Jrame comoving with the medium center of mass and ha ying four velocity 



itfi = (1, ). This is reflected in the second term of the RHS of (37 1. In this case, 



at T = for instance, e = £2/ + jJ.N and p = —Qf. Finally, when there is an exter- 
nal uniform magnetic field acting on the system, the additional symmetry breaking 
0(3) — > 0(2) takes place, and (t^ v ) get an anisotropy reflected in the appearance 
of the transverse metric structure rj^ v in (37 1. At T = we then have 



£ = £2f-fi 



d£2 



f 



2 : 



f 



-Qf + H- 



dn 



f 



B 2 



(38) 



(39) 



2 ' F " J ' " dB ' 2 
See Ref. ll66ll for detailed derivations of the formulas for the pressures and energy 
density in a magnetic field. 



8.2 MCFL Thermodynamic Potential 

Let us turn our attention now to densities large enough for the fermion system to 
be in the MCFL phase. Our ultimate goal is to find the EoS of this superconducting 
phase. To find the density and pressure of this phase, we first need, as seen from 
([3~8|-([39|, to obtain the contribution of the quarks to the thermodynamic potential. 
We can express the MCFL thermodynamic potential as the sum of the contributions 
coming from charged (Qc) an d neutral (Qn) quarks 11 131 

C2mcfl = Gc + Qn (40) 



with 



°c = ~Sz It 1 - ¥) r^ 3 e-^ +2 ^/A 2 [8 | e (c) | + 8 |g(c)| L (41) 

47r n=0 2 jQ 



Ci„ = ~[dp?e-?l>-\6\e<-<>->\+6\eW\\ 



and 



e^=±\(Jpl+2eBn-liy+Al 
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e (c) =±\l(Jpl + 2eBn + n) 2 + A 2 , (43) 



= ±yl(p-nY+Al ef = ±^{ P + nY + Al (44) 

being the dispersion relations of the charged (c) and neutral (0) quarks. In the above 
we used the notation 

k 2 _ Im _i_ . I \2 7 o A 2\2 



Kib = ~M±^ A2 +* A lY (45) 



The MCFL gaps A and As were introduced in Section 4. In the integrals (41 1 and 



(42 1, A -dependent smooth cutoffs for the NJL model are used. 

The effects of confinement can be incorporated by adding a bag constant S3 to 
Q.MCFL- Besides the bag constant and the quark contributions, the thermodynamic 
potential of the system also includes the pure Maxwell contribution, B 2 /2 l66ll . 
Hence, the thermodynamic potential of the MCFL phase is given by 

B 2 

n B = n MCFL + ,%+—, (46) 

The gaps A, and Ab have to be found from their respective gap equations 

d^MCFL _ „ dQ MCFL _ 

dA ' dA B ( ' 



8.3 EoS in a Magnetic Field 

The pressure and energy density of the MCFL phase are given by 

£mcfl = £2b-LI -j— , (48) 

Pmcfl = - q b, Pmcfl = -&b+B-^ (49) 

a B 

Note the splitting between parallel p^ MCFL (i.e. along the field) and transverse 
Pmcfl (i- e - perpendicular to the field) pressures due to the magnetic field. 



The magnetic field dependencies of the parallel and transverse pressures in (49 1 
were studied in Ref. [ 113], and are plotted in Fig. [2] Similarly to what occurs in the 
case of a magnetized uncoupled fermion system at finite density [66|, the transverse 



pressure in the MCFL phase increases with the field, while the parallel pressure de- 
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creases and reaches a zero value at field strength of order > 10 G for the density 
under consideration (fi = 500 MeV). We see from Fig.[2]that Qh and dCls/dB do 
not exhibit the Hass-van Alphen oscillations as happens with other physical quan- 
tities in the presence of a magnetic field l4Tll42ll46lllllllll21 . This is due to the 
high contribution of the pure Maxwell term in i2g and <9i2g/<9B, which makes the 
oscillations of the matter part negligible in comparison. 

1000-j- 
800 - 

m 
tg 600- 

5 

^ 400- 
200 - 
0- ■ 

1E17 1E18 1E19 

B (Gauss) 

Fig. 2 Parallel and perpendicular pressures as a function of the magnetic field intensity for fj. = 
500MeV and bag constant 38 = 5SMeV /fm 3 . 

The splitting between parallel and perpendicular pressures, shown in the vertical 
axis of Fig[3] grows with the magnetic field strength. Comparing the found splitting 
with the pressure of the (isotropic) CFL phase, we can address how important this 
effect is for the EoS. Notice that for 3 x 10 18 G the pressures splitting is ~ 10% of 
their isotropic value at zero field (i.e. the one corresponding to the CFL phase). 

1.0j 

0.8- 

? C6- 
m 

J, °- 4 - 

0.2- 
0.0- ■ 

1( 

B (Gauss) 

Fig. 3 Splitting of the parallel and perpendicular pressures, normalized to the zero value pressure 
(p(B = 0)), as a function of the magnetic field intensity for jl = 500 MeV and 3$ = 5&MeV/fm 3 . 
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In the graphical representation of the EoS in Fig.|4]the highly anisotropic behav- 
ior of the magnetized medium is explicitly shown. While the magnetic-field effect 
is significant for the e — p" relationship at B ~ 10 18 G, with a shift in the energy 
density with respect to the zero- field value of - 200 MeV/fm 3 for the same pres- 



sure, the field effect in the e - 
values. 



p relationship is smaller for the same range of field 




300 400 500 600 700 800 900 1000 



e (MeV/fin J ) 



300 400 500 600 700 800 900 1000 

e (MeV/fin 3 ) 



Fig. 4 Equation of state for MCFL matter considering parallel (right panel) and perpendicular (left 
panel) pressures for different values of B: zero field (solid line), 10 17 G (dashed line) and 5 x 10 18 G 
(dotted line). Note that the low value of B = 10 G is not distinguishable in the plots, being merged 
with the zero-field curve. The value of the bag constant was fixed to SS = SSMeV /fm 3 . 



The most important application of the EoS is to construct stellar models for com- 
pact stars composed of quark matter. This goal can be archived by using the rela- 
tivistic equations of stellar structure, that is, the well known Tolman-Oppenheimer- 
Volkoff (TOV) and mass continuity equations. 



dm , 
— = 4nre 
dr 



dP 
dr 



em 



4nr 3 P 



(50) 
(51) 



written in natural units, c = G = 1 . However, it is clear that this set of differential 
equations apply only to isotropic EoS for systems with spherical symmetry. 

If the magnetic field in the MCFL phase is high enough for the anisotropy in 
the pressure to be significant (i.e., expressed in terms of the pressure splitting to 
be (Ap/pcFL) ~ (B 2 /)J. 2 A 2 ) — the spherically symmetric TOV equations 

become inappropriate, because the deviations lead to significant differences with 
respect of realistic axi-symmetric models, yet to be constructed 11131 . 
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As we have stressed, an important characteristic of neutron stars is that they typ- 
ically possess very strong magnetic fields. Unveiling the interconnection between 
the star's magnetic field and the dense phases is important to understand the inter- 
play between QCD and neutron star phenomenology. As discussed above, in recent 
years much interesting work has been done on the properties of the different nuclear 
phases that can be reached in dense astrophysical objects in the presence of strong 
magnetic fields. An important new step in this context would be to consider the con- 
sequences for the star's phenomenology of the possible new phases. Although at this 
point we do not know yet the quantitative details of the potential consequences, no 
doubt exploring them will shed new light on the important question of how we can 
infer the presence of a color-superconducting core from astronomical observations, 
and whether such observations can distinguish among different color superconduct- 
ing phases. In what follow we discuss some of these related tasks. 



9.1 Low-Energy Physics 

The main challenge in determining the phases of matter inside a neutron star is to 
provide observables signatures of the presence of those phases. Currently there are 
many proposals to connect observations to the inner phases of the star. We want to 
discuss in general terms here, those connected to transport properties as conductiv- 
ities, viscosities, etc. As known, these transport properties are determined by the 
low-energy spectrum of the phase, that is, by the lowest-energy modes as Goldstone 
bosons and gapless quark excitations, which as already shown, can be affected by 
the presence of a sufficiently high magnetic field. Let's briefly mention some exam- 
ples of transport properties and how they could affect the star's observables. 

Viscosity. The viscosity of the interior of a star can be probed by observing how 
a rapidly spinning neutron star slows down. If the star slows down very quickly 
this indicates that it is unstable with respect to bulk flows (r-modes) that transfer 
the star angular momentum into gravitational radiation. But this can only occur if 
damping is sufficiently small. Based upon these arguments the possibility of pure 
CFL-quark-matter pulsars has been ruled out [51] since in the CFL phase viscous 
damping is negligible II 1031 . If the intermediate-density CS phase happens to have 
a large viscosity, it will not be restricted by r-modes arguments tough. 

As we showed in [56], in a three-flavor theory the spectrum of the NG bosons of 
the CS is affected by the restructuring of the gap produced by the magnetic field. 
As a consequence, instead of the 9 Goldstone bosons that exist in the CFL phase, 
in the MCFL only 5 remain. In contrast to the CFL case where several Goldstone 
bosons are Q-charged, in the MCFL all are Q-neutral. Therefore, the scattering rate 
of the low-energy bosons should be different in the magnetic background, and this 
will be reflected in turn in the transport properties of the star. By investigating trans- 
port properties as thermal conductivity and viscosity in the MCFL phase (or in the 
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extension of the MCFL phase when one takes into account the strange quark mass 
and neutrality effects) one could look for new observational effects that will allow 
us to distinguish between nuclear-core stars and quark-core stars. 

Thermal Conductivity. Given that neutrino emission rates and heat capacity gen- 
erally rise with density, neutron star cooling is likely preferentially sensitive to the 
properties of matter in the core. Investigations B1301I1311 on the impact of the ther- 
mal conductivity of dense quark matter on the star cooling process indicates that 
any CFL quark matter within the star will cool by conduction, not by neutrino emis- 
sion. In this direction, to investigate how the magnetic field can affect the medium 
thermal conductivity is of interest. 

Neutrino emission and detection . Neutrino emission is the dominant heat loss 
mechanism of the stars in their first million years. In ll84l 11231 the neutrino emis- 
sion from Nambu-Goldstone modes of the CFL phase has been investigated. These 
studies showed that the scattering of massless Goldstone modes, associated with 
the breaking of the baryon U(l)s symmetry, is not exponentially suppressed, and 
so, these modes dominate neutrino emission at late times. On the other hand, the 
time-of-arrival distribution of supernova neutrinos could be connected to possible 
phase transitions to and in quark matter |9l[30), but a detailed analysis of this sug- 
gestion requires a better understanding of both supernova itself and of the properties 
of quark matter at MeV temperatures. 

If the CS phase at intermediate densities results to be a variety of gapless phase, 
the gapless modes could play a significant role in the transport properties of the star. 
A recent study [9] has shown that even a relatively small region of gCFL matter 
in a star would dominate the heat capacity and the heat loss by neutrino emission. 
However, we already know that the gCFL is not stable, so it is unlikely that this 
phase will occur within the star. 

However, none of these studies have taken into account the presence of the in- 
medium magnetic field that penetrates the star's superconducting core. Neverthe- 
less, a total understanding of the transport mechanism of a compact star with a 
quark core will not be complete without considering the modification of the color- 
superconducting gap by the strong in-medium magnetic field, as well as by the mod- 
ification of the remaining Goldstone modes. 

This effect can be relevant for the low energy physics of a color superconducting 
star's core and hence for its transport properties. In particular, the cooling of a com- 
pact star is determined by the particles with the lowest energy; so a star with a core 
of quark matter and sufficiently large magnetic field can have a distinctive cooling 
process. This study is a pending task that is worth to be undertaken. 



9.2 Boosting Stellar Magnetic Fields via an Internal Mechanism 

The standard model 01371 11381 to explain the origin of the strong magnetic fields 
observed in the surface of magnetars is based on a magnetohydrodynamic dynamo 
mechanism that amplifies a seed magnetic field due to a rapidly rotating protoneu- 
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tron star. This model requires a spin period < 3ms. Nevertheless, this mechanism 
cannot explain all the observed features of the supernova remnants surrounding 
these objects |fl39l[T42l . 

As has been found recently, in color superconductors magnetic fields can be re- 
inforced l60l and even generated [63]. It is natural to expect that if a color super- 
conducting state exists in the core of neutron stars, it may have implications for the 
magnetic properties of such compact objects. At the moderately high densities that 
exist in the cores of neutron stars the most probable color superconducting state 
is not the CFL, but either an inhomogeneous phase or perhaps a strongly coupled 
2SC phase. In the 2SC phase, the Meissner mass of the charged gluons decreases 
with decreasing density to values which are close to zero. For such small charged 
gluon masses, a magnetic field does not need to be too large to induce a vortex state. 
Fields with values H > m 2 M will trigger the spontaneous generation of vortices of 
charged gluons, which in turn will enhance the existing magnetic field. Hence, CS 
could contribute to boosting the large magnetic fields observed in some stellar ob- 
jects as magnetars, without having to rely only on the quick spinning assumed in 
the standard model of magnetars that, as known, are associated with some of the 
conflict of this model with the observations. These induced gluon vortices could 
produce a magnetic field of the order of the Meissner mass scale, which implies a 
magnitude in the range ~ 10 16 — 10 17 G. Hence, the possibility of generating a mag- 
netic field of such a large magnitude in the core of a compact star without relying 
on a magnetohydrodynamics effect, can be an interesting alternative to address the 
main criticism 1 1391 11421 to the observational conundrum of the standard magne- 
tar's paradigm [83 92l II 091 fl 371 1 1 381 . On the other hand, to have a mechanism that 
associates the existence of high magnetic fields with CS at moderate densities can 
serve to single out the magnetars as the most probable astronomical objects for the 
realization of this high-dense state of matter. 



9.3 Stability of Magnetized Quark Stars 

It is now our goal to analyze the conditions for MCFL matter to become absolutely 
stable. This is done by comparing the energy density at zero pressure condition with 
that of the iron nucleus (~ 930 MeV). Depending on whether the energy density of 
the MCFL phase is higher or smaller than this value, the content of a magnetized 
strange quark could be or not made of MCFL matter. 

The stability criterion for MCFL matter can be derived in a simple way. Follow- 
ing Farhi and Jaffe's [45 1 approach, we can determine the maximum value of the 
bag constant that satisfies the stability condition at zero pressure for each magnetic 
field value. We call reader's attention that in all these derivations we work within 
a self-consistent approach, in which the solutions of the gap equations are substi- 
tuted in the pressures and energies of each phase before imposing the conditions of 
equilibrium and stability. 
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After imposing the zero pressure condition in the EoS for the MCFL phase, both 
the parallel and perpendicular pressures in the MCFL EoS need to vanish simulta- 
neously. Therefore, the two equilibrium conditions become 



Pmcfl 



Pmcfl — B 



&MCFL — 
d Q MCFL 

dB 



B- 



2 
~d~B 



0. 



B z = 



(52) 
(53) 



Where we are assuming that the bag constant depends on the magnetic field. It 
is not unnatural to expect that the applied magnetic field could modify the QCD 
vacuum, hence producing a field-dependent bag constant. One can readily verify 
that Eqs. (52 1-(53 i are equivalent to require p MCFL = and dp MCFL /dB = at the 
equilibrium point. 



Equation (53 1 can be rewritten as 



B=M- 



~d~B 



(54) 



where M = —dQ.MCFhl 'dB is the system magnetization. If we were to consider that 
the vacuum energy 3§ does not depend on the magnetic field, we would need 



M = B, 



(55) 



to ensure the equilibrium of the self-bound matter, a condition difficult to satisfy 
since it would imply that the medium response to the applied magnetic field (i.e. 
the medium magnetization M) is of the order of the applied field that produces it. 
Only if the MCFL matter were a ferromagnet this would be viable, but as known, 
the MCFL matter is on the contrary, an insulator. The other possibility for the equi- 
librium conditions p2]> and d53j to hold simultaneously is to have a field-dependent 



bag constant capable to yield nonzero vacuum magnetization Mq = — 



dB 



B. 




25 50 75 

M (MeV/fm 3 ) 



Fig. 5 Stability window for MCFL matter in the plane B vs. 2. The curve shown corresponds to 
the borderline value e/A = 930 MeV. 
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The following comment is in order. The fact that the bag constant needs to be 
field-dependent for self-bound stars in a strong magnetic field is a direct conse- 
quence of the lack of a compensating effect for the internal pressure produced by 
the magnetic field other than that applied by the vacuum (an exception could be of 
course if the paired quark matter would exhibit ferromagnetism). For gravitationally 
bound stars, on the other hand, the situation is different, since the own gravitational 
field can supply the pressure to compensate the one due to the field. For such sys- 
tems, keeping 38 constant in the EoS is in principle possible. Under this assumption 
we considered a fixed 38- value in Fig. [4] 

Taking into account that the matter energy density £ MCFL (i.e. the energy density 
that does not include the pure Maxwell contribution) divided by the baryon number 
is given by 

£ MCFL _ + & _ M d^MCFL 

n A n A n A djl 

We can write it under the zero parallel pressure condition (Qmcfl + 38 = —B/2) as 

+ (57 ) 
n A 2n A n A 

and the absolute stability condition becomes 

£ ^\, B =2^ + ^N<e Q (Fe 5 % (58) 
n A " 2n A n A 



Then, finding Hb as a function of B from (58 I, and substituting it back in (52 1, 
we can numerically solve 

38(B) = -n MCFL (n B ,B)-B 2 /2, (59) 

to determine the stability window in the plane B versus 38 for the MCFL matter 
to be absolute stable (Fig. [3). The inner region, which corresponds to smaller bag 
constants for each given B, is the absolutely stable region. 

Note that, contrary to Farhi and Jaffe RBI , we did not impose a minimum value 
for the bag constant because we have no clear indication from experiments of the 
possible behavior of this parameter when a magnetic field is applied to a system. 

In summary, our results indicate that a condition for the MCFL matter to be ab- 
solutely stable is a field-dependent bag constant that can give rise to a large vacuum 
magnetization at moderately strong fields (see Fig. |5j. Under these circumstances, 
increasing the magnetic field tends to destabilize the self -bound MCFL matter. This 
result differs from that found in [ 1 14] where it was used a CFL model at eB ^ 
with only one gap that was fixed by hand. 



Acknowledgements This work has been supported in part by DOE Nuclear Theory grant de- 
sc0002179 



32 



Efrain J. Ferrer and Vivian de la Incera 



References 



1. A. M. Abrahams and S. L. Shapiro, Astrophys. J. 374 (1991) 652. 

2. A. A. Abrikosov, Sov. Phys. JEPT 5 (1957) 1174. 

3. A. A. Abrikosov and L. P. Gorkov, Zh. Eksp. Teor. Fiz. 39 (1960) 1781 [Sov. Phys. -JETP 12 
(1961) 1243]. 

4. B. K. Agrawal, and S. K. Dhiman, Phys. Rev. D 79, 103006 (2009). 

5. M. Alford, K. Rajagopal and F. Wilczek, Phys. Lett. B 422 (1998) 247. 

6. M. Alford, K. Rajagopal and F. Wilczek, Nucl. Phys. B 537 (1999) 443. 

7. M. Alford, J. Berges and K. Rajagopal, Nucl. Phys. B 571 (2000) 269. 

8. M. Alford and S. Reddy, Phys. Rev. D 67 (2003) 074024. 

9. M. Alford, P. Jotwani, C. Kouvaris, J. Kundu, and K. Rajagopla, Phys.Rev. D71 (2005) 
114011. 

10. M. Alford, M. Brady, M. W. Paris, and S. Reddy, Astrophys. J. 629 (2005), 969. 

1 1 . M. Alford, et al., Nature (London) 445, E7 (2007). 

12. M. Alford, A. Schmitt, K. Rajagopal and T. Schcifer, Rev. Mod. Phys. 80 (2008) 1455 and 
references therein. 

13. M. Alford, Physics 3 (2010), 44. 

14. J. Ambjorn and P. Olesen, Nucl. Phys. B 315 (1989) 606. 

15. J. Ambjorn and P. Olesen, Phys. Lett. B 218 (1989) 67. 

16. P. W. Anderson, Phys. Rev. Lett. 3 (1959) 325. 

17. D. Bailin and A. Love, Phys. Rep. 107 (1984) 325. 

18. M. Baldo, et al., Phys. Lett. B 562 (2003) 153. 

19. V. Baluni, Phys. Rev. D 17 (1978) 2092. 

20. D. Bandyopadhyay, et al. Phys. Rev. D 58 (1998) 121301. 

21. B. C. Barrois, Nucl. Phys. B 129 (1977) 390. 

22. G. Baym, and S. Chin, Phys. Lett. B 62 (1976) 241. 

23. D. Blaschke, S. Fredriksson, H. Grigorian, and A. M. Oztas, Nucl. Phys. A 736 (2004) 203. 

24. M. Bocquet, S. Bonazzola, E. Gourgoulhon, and J. Novak, Astron. Astrophys. 301 (1995) 
757. 

25. A. Broderick, M. Prakash, and J. M. Lattimer, Astrophys. J. 537, 351 (2000). 

26. A. Broderick, M. Prakash, J. M. Lattimer, Astrophys. J. 537 (2000) 351. 

27. V. Canuto and H.-Y. Chiu, Phys. Rev,. 173 (1968) 1210. 

28. V. Canuto, and J. Ventura, Fund. Cosmic Phys. 2 (1977) 203. 

29. C. Y. Cardall, M. Prakash, and J. Ml Lattimer, Astrophys. J. 554 (2001) 332. 

30. G.W. Carter and S. Reddy, Phys.Rev. D62 (2000) 103002. 

31. R. Cassalbuoni, Z. Duan and F. Sannino, Phys. Rev. D 63 (2001) 114026. 

32. M. Chaichian, et al. Phys. Rev. Lett. 84 (2000) 5261. 

33. S. Chakrabarty, Phys. Rev. D 54 (1996) 1306. 

34. S. Chakrabarty, D. Bandyopadhyay, and S. Pal, Phys. Rev. Lett. 78 (1997) 2898. 

35. W. Chen, P.-Q. Zhang, and L.-G. Liu, Mod. Phys. Lett. A 22 (2007) 623. 

36. S. Coleman and D. Gross, Phys. Rev. Lett. 31 (1973) 851. 

37. J. C. Collins and M.J. Perry, Phys. Rev. lett. 30 (1975) 1353. 

38. P. Curie, J. Physique 3 (1894) 393. 

39. L. Dong and S. L. Shapiro, ApJ. 383 (1991) 745. 

40. A. Drago, A. Lavagno, and G. Pagliara, Phys. Rev. D 69 (2004) 057505. 

41. D. Ebert, K. G. Klimenko, M. A. Vdovichenko and A. S. Vshivtsev, Phys. Rev. D 61 (2000) 
025005. 

42. D. Ebert and K. G. Klimenko, Nucl. Phys. A 728 (2003) 203. 

43. E. Elizalde and E. J. Ferrer, and V. de la Incera, Ann. of Phys. 295 (2002) 33. 

44. E. Elizalde and E. J. Ferrer, and V. de la Incera, Phys. Rev. D 70 (2004) 043012. 

45. E. Farhi and R. L. Jaffe, Phys. Rev. D 30 (1984) 2379. 

46. B. Feng, E. J. Ferrer and V. de la Incera, Nucl. Phys. B 853 (201 1) 213. 

47. B. Feng, E. J. Ferrer and V. de la Incera, Phys. Lett. B 706 (201 1) 232. 



Magnetism in Dense Quark Matter 



33 



48. B. Feng, E. J. Ferrer and V. de la Incera, Phys. Rev. D 85 (2012) 103529. 

49. S. Ferrara and M. Porrati, Mod. Phys. Lett. A 8 (1993) 2497. 

50. E. J. Ferrer, V. de la Incera and A. E. Shabad, Phys. Lett. B 185 (1987) 407. 

51. E. J. Ferrer, V. de la Incera and A. E. Shabad, Nucl. Phys. B 309 (1988) 120. 

52. E. J. Ferrer, V. de la Incera, and A. E. Shabad, Annals Phys. 201 (1990) 51. 

53. E. I. Ferrer and V. de la Incera, Int. Jour, of Mod. Phys. A 11 (1996) 3875. 

54. E. J. Ferrer and V. de la Incera, Phys. Rev. D 58 (1998) 065008. 

55. E. J. Ferrer and V. de la Incera, Phys. Lett. B 481 (2000) 287. 

56. E. J. Ferrer, V. de la Incera and C. Manuel, Phys. Rev. lett. 95 (2005) 152002. 

57. E. J. Ferrer, V. de la Incera and C. Manuel, Nucl. Phys. B 747 (2006) 88. 

58. E. J. Ferrer, V. de la Incera and C. Manuel, PoS JHW 2005 (2006) 022. 

59. E. J. Ferrer, V. de la Incera and C. Manuel, J. Phys. A 39 (2006) 6349. 

60. E. J. Ferrer and V. de la Incera, Phys. Rev. Lett. 97 (2006) 122301. 

61. E. J. Ferrer and V. de la Incera, J. Phys. A: Math. Theor. 40 (2007) 6913. 

62. E. J. Ferrer and V. de la Incera, Phys. Rev. D 76 (2007) 04501 1. 

63. E. J. Ferrer and V. de la Incera, Phys. Rev. D 76 (2007) 1 14012. 

64. E. J. Ferrer and V. de la Incera, Phys. Rev. Lett. 102 (2009) 050402. 

65. E. J. Ferrer and V. de la Incera, Nucl. Phys. B 824 (2010) 217. 

66. E. J. Ferrer, V. de la Incera, J. P. Keith, I. Portillo and P. L. Springsteen, Phys. Rev. C 82 
(2010) 065802. 

67. For a review on the topic see M. Fiebig, J. Phys. D: Appl. Phys. 38 (2005) R123. 

68. E.S. Fradkin, Quantum Field Theory and Hydrodynamics, Proceedings of P. N. Lebedev 
Physical Institute, Vol. 29, 7 (1965, Moscow Nauka) (Engl. Transl., 1967 New York: Con- 
sultant Bureau). 

69. S. Frautschi, in: Proceedings of the Workshop on Hadronic Matter at Extreme Energy Density, 
N. Cabibbo, ed (Erice, Italy 1978). 

70. B. Freedman and L. McLerran, Phys. Rev. D 16 (1977) 1 166. 

71. B. Freedman and L. McLerran, Phys. Rev. D 17 (1978) 1 109. 

72. I. Fushiki, E. H. Gudmundsson, and C. J. Pethick, Astrophys. J. 342 (1989) 958. 

73. P. G. de Genness , Superconductivity of Metals and Alloys (W. A. Benjamin Inc, New York, 
1966). 

74. R. Gonzalez-Felipe, A. Perez-Martinez, H. Perez-Rojas, and M. Orsaria, Phys. Rev. C 77 
(2008) 015807. 

75. E. V. Gorbar, Phys. Rev. D 62 (2000) 014007. 

76. D. Gross and F. Wilczek, Phys. Rev. Lett. 30 (1973) 1343. 

77. V. P. Gusynin, V. A. Miransky and I. A. Shovkovy, Phys. Rev. Lett. 73 (1994) 3499. 

78. V. P. Gusynin, V. A. Miransky and I. A. Shovkovy, Nucl. Phys. B 563 (1999) 361. 

79. A.K. Hardings and D. Lai, Rep. Prog. Phys. 69 (2006) 2631. 

80. W. J. de Hass and P. M. van Alphen, Leiden Commun. A 212 (1930). 

81. W. J. de Hass and P. M. van Alphen, Proc. R. Acad. Sci. Amsterdam 33 (1930) 1106. 

82. J. Hufner, S.P. Klevansky, and P. Rehberg, Nucl. Phys. A 606 (1996) 260. 

83. Ibrahim et al, Astrophys. J. 609 (2004) L21. 

84. P. Jaikumar, M. Prakash, and T. Schafer, Phys.Rev. D 66 (2002) 063003. 

85. O. Kalashnikov and V. Klimov, Phys. Lett B 88 (1979) 328. 

86. W. Knzig, Ferroelectrics and Antiferroelectrics in F. Seitz, T. P. Das, D. Turnbull, and E. L. 
Hahn, Solid State Physics, 4 (Academic Press, p. 5, 1957). 

87. V.R. Khalilov, Phys. Rev. D 65 (2002) 056001. 

88. J. Kapusta, Nucl. Phys. B 148 (1979) 461. 

89. B. Keister, and C. Kisslinger, Phys. Lett. D 64 (1976) 1 17. 

90. K. G. Klimenko, Z. Phys. C 54 (1992) 323. 

91. K. G. Klimenko, Teor. Mat. Fiz. 90 (1992) 3. 

92. S. Kulkarni and D. Frail, Nature 365 (1993) 33. 

93. L. D. Landau, Z. Phys. 64 (1930) 629. 

94. L. D. Landau and E. M. Lifshitz, The Classical Theory of Fields Elsevier Butterworth- 
Heinemann (Four Edition, 1975) MA, Ch. 11 (pp. 290-293). 



34 



Efrain J. Ferrer and Vivian de la Incera 



95. L. D. Landau and E. M. Lifshitz, "Electrodynamics of Continuous Media" 1960 (Oxford: 
Pergamon). 

96. D. S. Lee, C. N. Leung and Y. J. Ng. Phys. Rev. D 55 (1997) 6504. 

97. C. N. Leung, Y. J. Ng and A. W. Ackley, Phys. Rev. D 54 (1996) 4181. 

98. A. D. Linde, Phys. Lett. B 86 (1979) 39. 

99. M. Lines and A. Glass, Principles and applications of ferroelectrics and related materials, 
(Clarendon Press, Oxford, 1979). 

100. D.F. Litim and C. Manuel, Phys. Rev. D 64 (2001) 094013. 

101. G. Lugones and J. E. Horvath, Phys. Rev. D 66 (2002) 074017. 

102. G. Lugones, and I. Bombaci, Phys. Rev. D 72 (2005) 065021. 

103. C. Manuel, A. Dobado and F.J. LLanes-Estrada, JHEP 0509 (2005) 076. 

104. D. P. Menezes, et al., Phys. Rev. C 79 (2009) 035807. 

105. D. P. Menezes, et al., Phys. Rev. C 80 (2009) 065805. 

106. P. Morley, and M. Kisslinger, Phys. Rep. C 36 (1979) 3. 

107. N. F. Mott, Rev. Mod. Phys. 40, 677 (1968). 

108. N. F. Mott, Metal-Insulator Transitions (Taylor and Francis, London, 1974). 

109. T. Murakami et al, Nature 368 (1994) 127. 

110. N. K. Nielsen and P. Olesen, Nucl. Phys. B 144 (1978) 376. 

111. J. L. Noronha and I. A. Shovkovy, Phys. Rev. D 76 (2007) 105030. 

112. K. Fukushima and H. J. Warringa, Phys. Rev. Lett. 100 (2008) 03200. 

113. L. Paulucci, E. J. Ferrer, V. de la Incera, and J. E. Horvath, Phys. Rev. D 83 (201 1) 043009. 

1 14. A. Perez Martinez, R. Gonzalez Felipe and D. Manreza Paret, Int. J. Mod. Phys. E 20 (201 1) 
84. 

115. A. Perez-Martinez, H. Perez-Rojas, and H. J. Mosquera-Cuesta, Eur. Phys. J. C29 (2003) 
111. 

116. H. Perez Rojas and A. E. Shabad, Ann. Phys. (N. Y.) 121 (1979) 432. 

117. G. Pagliara, and J. Schaffner-Bielich, Phys. Rev. D 77 (2008) 063004. 

118. D. Page, and S. Reddy, Ann. Rev. Nucl. Part. Sci. 56 (2006) 327. 

119. H. Politzer, Phys. Rev. Lett. 30 (1973) 1346. 

120. A. Rabhi, et al., J. Phys. G 36 (2009) 115204. 

121. A. Rabhi, C. Providencia, and J. Da Providencia, J. Phys. G 35 (2008) 125201. 

122. R. Rapp, T. Schafer, E. V. Shuryak, and M. Velkovsky, Phys. Rev. Lett. 81 (1998) 53. 

123. S. Reddy, M. Sadzikowski, and M. Tachibana, Nucl. Phys. A 714 (2003) 337. 

124. Rgnvaldsson et al., Astrophys. J. 416 (1993) 276. 

125. W. C. Rontgen, Ann. Phys. 35 (1888) 264. 

126. S. B. Ruester, and D. H. Rischke, Phys. Rev. D 69 (2004) 04501 1. 

127. A. Schmitt, Q. Wang, and D. H. Rischke, Phys. Rev. D 69 (2004) 094017. 

128. S. L. Shapiro and S. A. Teukolsky, Black Holes, White Duarfs, and Neutron Stars (Wiley, 
New York, 1983). 

129. B. K. Sharma, P. K. Panda, and S. K. Patra, Phys. Rev. C 75 (2007) 035808. 

130. I. Shovkovy and P.J. Ellis, Phys. Rev. C 66 (2002) 015802. 

131. I. Shovkovy and P.J. Ellis, Phys. Rev. C 67 (2003) 048801. 

132. V. V. Skalozub, Sov. J. Nucl. Phys. 23 (1978) 113. 

133. V. V. Skalozub, Sov. J. Nucl. Phys. 43 (1986) 665. 

134. D. T. Son, Phys. Rev. D 59 (1999) 094019. 

135. D. T. Son and M. A. Stephanov.Phys. Rev. D 61 (2000) 074012. 

136. I.-S. Suh, and G.J. Mathews, Astrophys. J. 546 (2001) 1126. 

137. C. Thompson and R. C. Duncan, Astrophys. J. 392 (1992) L9. 

138. C. Thompson and R. C. Duncan, Astrophys. J. 473 (1996) 322. 

139. J. Vink and L. Kuiper, Mon. Not. Roy. Astron. Soc. Lett. 370 (2006) L14. 

140. FX. Wei, et al. J. Phys. G 32 (2006) 47. 

141. H. A. Wilson, Phi l. Trans. R . Soc. A 204 (1905) 129. 

142. R-X Xu, [astro-ph/0611608"l 

143. F. Yang, H. Shen, Phys. Rev. C 79 (2009) 025803. 



